Mapping Integer Order Neumann Functions 
g: To Real Orders 
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Abstract 

In a recent paper we unified Bessel functions of different orders .Here we extend 

f- >, | the unification to other linairely independant solutions to Bessel equation ,Neu- 

r , "x. ■ mann's and Hankel's functions 
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Neumann's functions N p (z) and Bessel functions J p (z) with p integer ,are the two 
linairely independant solutions to the Bessel equation |l| . 

z 2 y"(z) + zy\z) + (z 2 - p 2 )y(z) = (1) 

When p is real it is known that J p (z) and J- P (z) are also two linairely independant 

solutions of the above equation so consequently we can rewrite N p (z) as a combination of 

Bessel functions of positive and negative orders and the combination is as follows 

NJz) = J M cos ^ ~ J -^ ( 2 ) 

sinjm 

Bessel and Neumann's functions are the "analog" respectively of the cosine and sine ,while 
Hankel functions H^p and H^ yet another linear combinations of J p (z) and J~ p (z) are 
the "analog" of e ipz and e~ ipz . Hankel functions are related to N p (z) and J p (z) as. 
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N p (z) = H » {Z) 2i H * {Z) (3) 

where H^- 1 ' and H ( 2 > are complex conjugate to each other. 

In a recent paper || we showed, for the first time , that Bessel functions of real orders 
J n +x < A < 1 and Bessel functions J n n G Z of integer orders are no longer disassociated 
objects , but one is the deformedQ version of the other and as a consequence ,they are 
related to each other through the unifying formula 

■ 7 " +aU) ez P (-AE)^# (4) 
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At this point it seems natural to extend this formula which applies to J n , to the second 
linairely independant solution iV n .The apparently straigthforward way to do would be to 
use formula [| |5| and then apply the mapping operator e ±A ^ , but unfortunately our 
formula applies to the reduced Bessel function J n (z)/z n rather than directly to J n (z) and 
this make our formula of no practical use for our purpose. In the following we will propose 
a nice way to relate N n+ \ to N n . 

lr The deformation mechanism is described in B 



Let us recall first that to Bessel functions one can associate polynomials O n (z). These 
are Neumann polynomials .The relationship between them is performed through the 
known formula. 



J2Uz)A n (t) (6) 
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with O n (z) = A n + (-l)M_ n 

The first step towards the solution of our problem is to show that if , to Bessel functions 
one can associate Neumann polynomials O n (z), one can equally well associate to reduced 
Bessel functions something and these associated objects are precisely Neumann function 
N n (z) . We are not sure if a formula which relates reduced Bessel functions to Neumann 
functions , similar to the above one |6], already exists in the literature ,so we prefer to work 
it out first. In doing this we largely follow a method which dates back to Sonine [[| and 
adapts it to the present need. Let us first state the result. Let fi(r) be an an arbitrary 
function of r and if ^(t) = x ,let r = U(x) so that 15 is the inverse function to Q .Denote 
the function and the function we associate to it respectively Z n (z) and A n (t) and define 
them through their integral representations as follows 

ZJz) = -, i e~* 2n+ i </; 
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A n {t) = / e l x ~ a {U{x)) n dx 
Jo 

then it is not that difficult to prove the existence of a formula relating both functions 
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(7) 



To prove the result [7| , suppose for any given positive value of x that | r |>| U(x) \ on a 
closed curve C surrounding the origin and the point z and that | r |<| U(x) \ on a closed 
c surrounding the origin but not enclosing the point z .Then compute the series 

£ Z n A n = — Ef / exp(-z 2 n(r)+t 2 x + -(T-U(x))^ r dxdT 

n=-oo ^ nl n=0 JeJo IT 

I +°° r roo \ T n 

+ V & / exp(-z 2 tt(T) +t 2 x + -(t-U(x)) , r - dx dr 

2m^ JcJo n u 2 V K JJ U^x)^ 1 

1 r°° r rexp(-z 2 VL(T)+t 2 x + \{r-U(x) 
2~jxik Vc~ tc r-U(x) T 

1 r 00 r exp(-z 2 Q(r) + t 2 x + \{r - U(x) 
2iri Jo Jz t — U(x) 
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In going to the third line,we perform the summations ^ J2t=o \~ t) = T -u(x) ( r on ^ e 
path C )and ^ J2n=o \uh)) = ~ T -u(x) ( r on ^ e P a ^ n c ) which are convergent according 
to the conditions above .For the case of interest we have Z n (z) = - J 4r^ , fi(r) = ^- and 
U(x) = w- -In this case the function associated to the reduced Bessel function has the 
integral form . 
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The integral defining the function A n (t) is easily identified to the K n (t) function 
which is the Hankel function of imaginary argument K v (t) = ^e % ^ v% H. _ v {it) . 

i , t 2 x dx K-n + i(it) K_ n+ i(—it) 
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2 Jo x n (it)- n+1 (-it)- n+1 

with | arg t |< 7r/2 and Re t 2 > 

Equation [7| then relates reduced Bessel funnctions to Hankel functions which are as- 
sociated to them 



2 l +5? JJz 



-lH-(l) 



7T £ 2 — Z 2 _ 2 



E^^^n-iW- (9) 
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In writing this formula we have used the property of the Hankel functions that H- n = 
(-l) n H n . 

At this point if we restrict ourselves to t and z reals and take the complex conjugate 
of the above formula we get another formula 

-g^r-^iW- (io) 
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where, this time, we use the property Hj^(z)* = Hf®(z*) . Adding || to [II] and using the 
definition of Neumann function in term of Hankel functions |3| we get the desired result 

1 +00 T / \ 

E^^-^n-lW. (11) 



7Tt 2 ~Z 2 ^ 

The particular form of the left hand side of this equation (a happy event) and the presence 
of reduced Bessel functions on the right hand side of the equation allow a straigthforward 
and quite elegant application to Neumann 's function of the unifying formula f| ( valid for 
reduced Bessel functions ) .In fact we can perform a succession of valid operations . 

^(-ADj^ = E^HD^^-iW (12) 



2 ^(-aE)A = E^^^iWt) 



T t 2 - z 2 ~ z n+x 



l W ^ = E^artAD*"- 1 ^*) 
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Where the capital sigma operator with the t variable acquires an extra factor (— l) m . 

E=E(-ir^ as) 

t meZ/O 

This extra factor comes from the identity. 
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dz 2 t 2 -z 2 dt 2 t 2 -z 2 



(14) 



Note that in the last line of [12| we identified the result of the action of the operator 
exp(X J2t) on Neumann functions of integer orders with Neumann function of orders n + X 
as for reduced Bessel functions. That is we put 

t n+x N n+x (t) = exp(\ J2) t n N n (t) (15) 



It remains to show [15|, that is ,the result of the mapping operation by the operator 
e x J2 is indeed Neumann functions of real orders n + X and not something else. To see this 
let us proceed as follows 
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In the first line we use the fact that the sum is independant of A and hence change the 
parameter A —* X — m .In the second line we absorb the index m as we sum over all m e Z 
.From |T^ we infer the symmetry property of the A n \ 

A„ iA = A„ +mjA _ m V A, m (17) 

This symmetry means that A n \ only depends on the combination n+X that is A n \ = A n+ \ 
.This remark together with the fact that A n+X \x=o— t n N n (t) leads to the identification 



15 as we anticipated in the last line of 12 



Formula |T5| is the main result of the paper .It shows that Neumann functions are unified 
irrespective of their orders being integer or real .They are unified through the same formula 
as reduced Bessel functions z n J n (z) in |5|. 

Note that we could have already applied the procedure to Hankel function's first by 
starting from the identities ||] and . [To| . In this case we could have obtained a similar 
unifying formulas ( i=l,2 ). 

t n+X H { :Ut)=exp(Xj:)t n H^(t) (18) 



Let us note that our formula ^ for the evolution of Hankel functions could be used to infer 
the evolution equation for reduced Bessel function z n J n by using the defining equation |3] 
of J ,s in terms of H' s . By simple inspection one can see that the result we get is indeed 
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